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On Legendre's Function P„(#), when n is great and has any 

Value. 

By Lord Kayleigh, O.M., F.E.S. 
(Received April 27, 1916.) 

As is well known, an approximate formula for Legendre's function T n (8) } 
when n is very large, was given by Laplace. The subject has been treated 
with great generality by Hobson,* who has developed the complete series 
proceeding by descending powers of n, not only for P n but also for the 
" associated functions." The generality aimed at by Hobson requires the 
use of advanced mathematical methods. I have thought that a simpler 
derivation, sufficient for practical purposes and more within the reach of 
physicists with a smaller mathematical equipment, may be useful. It had, 
indeed/ been worked out independently. 

The series, of which Laplace's expression constitutes the first term, is 
arithmetically useful only when nO is at least moderately large. On the 
other hand, when 6 is small, V n tends to identify itself with the Bessel's 
function Jo(nd) } as was first remarked by Mehler. A further development 
of this approximation is here proposed. Finally, a comparison of the results 
of the two methods of approximation with the numbers calculated by 
A. Lodge for n = 20f is exhibited. 

The differential equation satisfied by Legendre's function P n is 

ff 7/ flit * 

d0* + ' GOt0 a§ + n ( n + 1 ) u = °- (!) 

If we assume u = v (sin 0)-*, and write m for n + $, we have 

+ m 2 v = — r-— . (2) 

If we take out a further factor, e im9 } writing. 

u = v sin"* 6 = w e im6 sin~s #, (3) 

of which ultimately only the real part is to be retained, we find 

dho , . dw . w A 

dff 2 d6 4sm 2 6> v { 

*"Ona Type of Spherical Harmonics of Unrestricted Degree, Order, and Argu- 
ment," 'Phil. Trans./ A, vol. 187 (1896). 

+ " On the Acoustic Shadow of a Sphere," < Phil. Trans.,' A, vol. 203 (1904) ; * Scientific 
Papers/ vol. 5, p. 163. 
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We next change the independent variable to z, equal to cot 9, thus 
obtaining 

Tz-2^\} 1+Z) dz> +2 *fc + ZS m . (0) 

From this equation we can approximate to the desired solution, treating m as 
a large quantity and supposing that to = 1 when z — 0, or 9 = |7r. 
The second approximation gives 

dw i -, . ^ 

whence w = 1 — 



dz 8 m 8m* 

After two more steps we find 

W ^V8m 128m 8 / 128m 2 1024m 3 ' {)> 



Thus in realised form a solution of (1) is 

9 cot 2 ff 
128m s 



% = C sin * 



1 — -^-^ — - ) cos (me/ + 7) 



, fcot 9 9 cot 9 75 cot 3 01 . . , A , x ~| /H7X 

Lorn L28m 3 1024 m 6 J J 

and this may be identified with P n provided that the constants C, 7, can be 
so chosen that u and dujdO have the correct values when 9 = -|7r. For this 
value of 9 we must have 

P* (i tt) = C cos (I- mTr + 7), (8) 

(dP./d%= c(-m- ^ + I J^)sinam7r + 7 ). (9) 

We may express {clV n jcW\„ by means of Pk+iQtt). In general 

sin 2 -^ = (» + 1) (cos ^P.-Ph+i), 

Cv COS 1/ 

so that when 9 = ^7r, 

<£P w /d0 = -d? n /d cos = (% + l)P»+i. (10) 

When % is ^<m, (d'P n /d9)i 7r vanishes, and, C being still undetermined, we 
may take to satisfy (9), 7 = — \ir\ and then from (8) 

so that 

p _ 1 . 3 . 5...(%— 1) 

2.4.6... n 

Here n is even, say 2r, and it is supposed to be great. Thus 

r 1.2.3...(2r-l) 2r (2r)i . 

2 2 .4 2 .6 2 (2r) 2 2 2p (r'!) 2 ' 
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and when r is great, 



Thus 



c = -T7— r Ex P- ( hi: :- 0^-^ + 



v/(7rr) r \24r 360. 8r 3 6r 360r 3 , 



v/(7rr)\ 8r 128 r* 1024r 3 / 



\.7m/ \ 4% 32 % 2 128 ?i 3 
When % is even and with this value of C, 



P M = C sin 2 51 



, /eot# 9eot# 75cot 3 #\ . f/ , 1XZ3 -, ^ /10N 

When w is odd, the same value of 7, viz. —£77-, secures the required 
evanescence in (8), and we may conjecture that the same value of C will also 
serve. Laplace* indeed was dbntent to determine 7 from the case of n odd 
and C from the case of n even. I suppose it was this procedure that 
Todhunterf regarded as unsatisfactory. At any rate there is no difficulty in 
verifying that (9) is satisfied by the same value of C. From that equation 
and (10), 

C = ( - 1 W+V ( n + 1 ) p n+i( i i rr) 
v ; m + l/8m-9/128m 3 ' 

and 

.L.O.O... Tt* 



(-!)»(»•") P» +1 (iir) 




2. 4. 6... 0+1) 
2 1 ^_ 1 + l^ 



.(^4-l)7rJ L 4(^+1) 320i+l) 2 128(^+1) 3 . 

Here, as throughout, m = ^4-J, and when we expand these expressions in 
descending powers of n we recover (11). Equations (11) and (12) are thus 
applicable to odd as well as to even values of n. 

But whether n be even or odd, (12) fails when 6 is so small that n6 is not 
moderately large. For this case our original equation (1) takes approxi- 
mately the form 

dP+m + * tt = > (13) 

* * M6c. Cel., J Supplement au V e volume. 
t ' Functions of Laplace, etc.,' p. 71. 
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where a 2 is written for n(n + l); and of this the solution is 

n = J {ad). (14) 

It is evident that the Bessel's function of the second kind, infinite when 
6 = 0, does not enter, and that no constant multiplier is required, since u is 
to be unity when = 0. For a second approximation we replace (13) by 

d 2 u . Idu , 9 du/l cos#\ flcifo aO T , , m 
<i# 2 d0 c£# \tf sm 0/ 3 <i# 3 k 



or, if «# = z } 



cl 2 u , ldu 
7l 



2 + zfa +U = —* J °'(*)- 



Z 

3 a* 



(15) 



In order to solve (15) we assume as usual 

'u — v.J Q (z). 



This substitution gives 



d?vdv (2 Jq' ] 1\ _s__.y 
^ 2 3 ^ \ Jo «/ 3 a? Jo 



(16) 



(17) 



a linear equation of the first order in dvfdz. In this 

j ("j7 + j)^ = lo g(^Jo 2 ); 



so that 



dv _ A 1 

c^ 2J0 2 3a 2 a;Jo 



£ 2 J Jo'dz . 



«JL JLp v-/ -A. V/ 



J^JoJo'cfe = ^ 2 Jo 2 - J Jo 2 ^ = KJo 2 -ia 2 (Jo a + Jo' 2 ) = 



1 ,2T '2 



X XI US 



$y A 



T 'i 

on 



3 tlO 



ete zJo 2 6 a 3 J 2 ' 



(18) 



which has now to be integrated again. 



J /2 £ dz 

17~~ 



— * I zJo'd[-j-) 
\<J o/ 



^Jo . r ^Jo +Jq ^ 



Jo 






Jo 



sJV 

' Jn 



zdz = 



« Jo' s 2 



regard beins 



jLXlLlo 



paid to the differential equation satisfied by J . 
v = B-f A 



'dz 1 

Jjn~ 2 " + "l2^^ +2 ^°/ Jo)> 



and 



^ =■ B Jo + A Jo 



cfe 1 

"TTiTT^ 5(3 Jo + 2 2 Jq ). 

2 Jo 2 12 a 2, v y 



(19) 



(20) 



when n is great and 6 has any Value. 
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For the present purpose A = 0, B = 1 ; so that for P w , identified with % 



we get 



P n (0) = J (z) + -f~ {z*J (z) + 2zJ ' (z)} , 



12 a? 



(21) 



; = a$ y a 2 = n (n + 1). 
Ji, are thoroughly tabulated.* 



in which 

The functions J , Jo' = 

The Table annexed shows in the second column P20 calculated from (21) 
for values of ranging from 0° to 35°. The third column gives the results 
from (11), (12), beginning with = 10°. In the fourth column are the 
values of P« calculated directly by A. Lodge. It will be seen that for 
= 15° and 20° the discrepancies are small in the fifth place of decimals. 
For smaller values of 0, the formula involving the BessePs functions gives 
the best results, and for larger values of the extended form of Laplace's 
expression. "When exceeds about 35° the latter formula gives P20 correct 
to six places. For n greater than 20 the combined use of the two methods 
would of course allow a still closer approximation. 



Table for P 2 q. 



e. 


Formula (21). 


From (11) and (12). 


Calculated by Lodge. 






1 -oooooo 




1 -oooooo 


5 


'346521 


— 


-346521 


10 


-0-390581 


-0-390420 


-0*390588 


15 


-0-052776 


-0-052753 


-0-052772 


20 


+ -300174 


+ 0-300191 


+ 0-300203 


25 


-0-078051 


-0-078085 


-0-078085 


30 


-0-216914 


-0-216997 


-0-216999 


35 


+ 0-3 55472 


+ -155635 


+ 0-155636 


40 





+ 0-127328 


+ -127328 


45 


— — 


+ -193065 


+ -193065 



* See Gray and Mathew's * BessePs Functions. 1 
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